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Abstract 

The connection between noncentrosymmetric materials’ structure, electronic structure, and bulk 
photovoltaic performance remains not well understood. In particular, it is still unclear which pho¬ 
tovoltaic (PV) mechanism(s) are relevant for the recently demonstrated visible-light ferroelectric 
photovoltaic (K,Ba)(Ni,Nb) 03 _( 5 . In this paper, we study the bulk photovoltaic effect (BPVE) of 
(K,Ba)(Ni,Nb) 03_(5 and KNb 03 by calculating the shift current from first principles. The effects of 
structural phase, lattice distortion, oxygen vacancies, cation arrangement, composition, and strain 
on BPVE are systematically studied. We find that (K,Ba)(Ni,Nb) 03_(5 has a comparable BPVE 
with that of the broadly explored BiFe 03 , but for a much lower photon energy. In particular, the 
Glass coefficient of (K,Ba)(Ni,Nb )05 in a simply layered structure can be as large as 12 times that 
of BiFe 03 . Furthermore, the nature of the wavefunctions dictates the eventual shift current yield, 
which can be significantly affected and engineered by changing the O vacancy location, cation 
arrangement, and strain. This is not only helpful for understanding other PV mechanisms that 
relate to the motion of the photocurrent carriers, but also provides guidelines for the design and 
optimization of PV converters. 
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I. INTRODUCTION 


As the world power consumption and carbon emissions continue to increase, solar energy 
has drawn even more attention because it is clean, abundant, and sustainable, and thus is 
widely seen as a long-term substitute for traditional fossil fuels [I]. ElRcient solar energy 
conversion relies primarily on semiconducting photoabsorbers with a low band gap (Eg), 
allowing absorption of most of the solar light reaching the earth. With light absorption, 
electrons are excited to the material’s conduction band (CB) for electricity generation or 
catalysis [2H1!. However, the photo-excited electrons may also recombine radiatively or non- 
radiatively with the created holes, reducing the power conversion efficiency. In conventional 
solar cells, this recombination rate is minimized by an externally engineered asymmetry, i.e., 
electrons and holes are separated by the electric field in a p — n junction or other interface. 
This not only complicates the device fabrication, but also imposes the Shockley—Queisser 
limit on the power conversion efficiency of this type of device [5]. Ferroelectric (FE) solar 
converters, lacking inversion symmetry due to intrinsic spontaneous polarization, can sep¬ 
arate photo-excited charges by the depolarization field or by the bulk photovoltaic effect 
(BPVE) [HHH]- In the BPVE, a spontaneous direct short-circuit photocurrent is generated 
when electrons are continuously excited to a quasiparticle coherent state that has an intrinsic 
momentum, avoiding the need for an interface to separate charge. In particular, the BPVE 
is able to generate an above-band-gap photovoltage [9] , potentially enabling a higher power 
conversion efficiency than the Shockley—Queisser limit. 

However, most conventional ferroelectric materials have wide band gaps [Eg >2.7 eV for 
BiFeOs, Eg >3.5 eV for Pb(Zri/ 2 Tii/ 2 ) 03 ], limiting their ability to absorb the visible light 
that makes up the biggest fraction of the solar irradiance. Thus, an enormous amount of 
efforts has been focused on the design and optimization of FE materials in order to reach 
a lower band gap HUDM. Among them, the study and improvement of ferroelectric ox¬ 
ides are particularly important, as these materials can be integrated into the conventional 
electronics [2BHH1]- In particular, recently a ferroelectric perovskite (K,Ba)(Ni,Nb) 03_5 
(KBNNO) has been successfully synthesized and demonstrated to have a near-optimal band 
gap (1.39 eV), exhibiting good photovoltaic (PV) performance [35]. While this has sub¬ 
stantially advanced the area of ferroelectric photovoltaics, there remain open questions. For 
example, what is the underlying PV mechanism in this material? Furthermore, this material 
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is also able to exhibit BPVE, as previous time dependent perturbation theory analysis has 
shown that BPVE, in principle, can arise in any polar material through the “shift current” 
mechanism [HSIIHT!. Actually this becomes even more fascinating when taking into account 
that the parent KNbOa is an interesting nonlinear optical (NLO) material with high NLO 
coefhcients |38l|39]. In this paper, we study the BPVE and its correlation to structural and 
electronic properties in KBNNO and KNbOa from hrst principles. The connection between 
the photocurrent and electronic structure elucidated here is not only useful for understand¬ 
ing other PV mechanisms in KBNNO in the sense that all these PV mechanisms relate to 
the light absorption and the motion of the photo-excited carriers and thus to the electronic 
properties, but also can be generalized to other similar materials that have great tunability 
of orbital character near the band gap. 

II. COMPUTATIONAL DETAILS 

Density functional theory (DPT) calculations were performed with the local density ap¬ 
proximation (LDA) functional, as implemented in the QUANTUM-ESPRESSO code [101 - 
112]. Norm-conserving, optimized nonlocal pseudopotentials were used to represent all ele¬ 
ments [H]. The DPT-I-I7 method was used to improve the description of d-orbital electrons 
by better accounting for the correlation effect, with Hubbard U parameterized by the linear- 
response approach m. The calculated Hubbard U values are 3.6 eV for Nb in KNbOa, 
and 3.97 eV and 9.90 eV for Nb and Ni in KBNNO. The shift current was calculated with 
a previously developed first-principles approach based on time-dependent perturbation the¬ 
ory, which yields good agreement with experiment for the prototypical ferroelectric oxides, 
such as BiEeOs [HTIIIS]. To calculate the shift current, a self-consistent calculation was first 
done by the LDA-I-U approach for the structures fully relaxed by LDA, followed by non-self- 
consistent calculations with much finer k grids. We used the LDA relaxed structures for shift 
current calculations, as LDA was shown to describe well the KNbOs structural properties, 
with only 0.3% underestimation of the tetragonality c/a [23|- The Monkhorst-Pack A:-mesh 
method [IH] was used to sample the Brillouin zone. To converge the shift current, the k 
grid must be sufficiently dense, e.g., a 40x40x40 k grid was used for a typical ABO 3 unit 
cell. In order to overcome the self-interaction error of the standard DPT method, we also 
used the HSE06 hybrid functional to calculate the band gaps of some compositions [17| . The 
HSE06 hybrid functional improves the band gap description by including a proportion of the 
exact exchange interaction, while the correlation part remains the same as in the standard 
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DFT method; it only serves as a corroboration of the band gap here. We also adopted the 
GGA+U method to calculate the shift current of KNbOs and some KBNNO structures, and 
confirmed that the shift current spectral features are essentially similar to those obtained by 
the LDA+G method (the change of the shift current magnitude is tiny), except for slightly 
larger onset photon energies in the GGA+G case. 

III. RESULTS AND DISCUSSIONS 
A. The BPVE of KNbOa 

Shift current, as a dominant mechanism for the BPVE, is a second-order nonlinear optical 
effect with the photocurrent quadratic in the electric field {Jg = arsqErEg) |37j. The Glass 
coefficient 


G 


rrq 
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describes the current response in a thick sample and includes the light attenuation effect 
due to the absorption coefficient 

We start by calculating the shift current of the parent material KNbOa. KNbOa is 
a typical ABO 3 perovskite ferroelectric oxide that occurs in various different phases. At 
temperatures above 691 K, KNbOs is in a paraelectric cubic phase with space group Pm3m. 
As the temperature decreases below 691 K, it first undergoes a phase transition into a 
tetragonal phase (space group P4mm), and then into an orthorhombic phase at 498 K (space 
group Amm 2 ) and an even more distorted rhombohedral phase with space group R 3 m at 



FIG. 1. (Color online) The (a) shift current susceptibility and (b) Glass coefficient of various phases 
of KNbOs. 
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TABLE I. The HSE06 band gap of various KNbOs phases (cubic, tetragonal, orthorhombic, 
and rhombohedral) and two different cation arrangements (1x1x2 and rocksalt) of the 
(K,Ba)(Ni,Nb)05 (KBNNO) solid solutions. KBNNO has much smaller band gaps than KNbOa. 
All structures are fully relaxed with LDA. 


System 

KNbOs 

(K,Ba)(Ni,Nb)05 

Cubic Tetragonal Orthorhombic Rhombohedral 

1x1x2 rocksalt 


pHSE 

9 

2.65 2.71 3.10 3.33 

1.28 2.16 


263 K. Concurrently with the structural transitions from cubic to tetragonal, orthorhombic, 
and rhombohedral phases, the NbOe octahedra become more severely distorted through 
octahedral tilting and rotation, and the Nb ions move away from the center of the Oe cages. 
As a typical ABO 3 perovskite semiconductor/insulator, the valence band maximum (VBM) 
of KNbOs arises mainly from the O 2 p orbitals, while the conduction band minimum (CBM) 
is predominantly composed of Nb 4d orbitals [IHl HH] • Also, the band gap of KNbOs increases 
significantly when the structure changes from cubic to the more distorted rhombohedral 
phase, because of the correlation between the NbOe octahedral distortion and Nb ions off- 
center displacements and the electronic structure |22]. Correspondingly, experiment finds 
a wide range of band gap values, from 3.3 eV for the cubic phase [50] to 4.4 eV for the 
tetragonal phase Our previous HSE06 calculation also shows that the band gap of the 
rhombohedral KNbOs is ~0.6 eV larger than that of its tetragonal counterpart (Table|I]) [22|- 

Figure [T] shows the calculated shift current susceptibilities and Glass coefficients of KNbOs 
in its ferroelectric tetragonal, orthorhombic, and rhombohedral phases. For convenience, we 
only shows the largest tensor element xxZ, where the upper case letter represents the shift 
current direction. Clearly, the magnitude of the maximum shift current susceptibility for 
the room-temperature orthorhombic phase with respect to the photon energies between 
Eg and 4-1.0 eV (~15 xlO ^ V is more than twice that of the more broadly studied 
BiFeOs (!5i6xl0“^ V“^), although their band-edge shift current responses are on a par with 
each other (?alxl0“^ V“^). This indicates that KNbOs is more promising than BiFeOs 
for photovoltaic applications with high photon energies (ultraviolet light), but it is not as 
good as BiFeOs for light at the visible-UV edge (3.0-3.2 eV), since BiFeOs has a much 
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lower band gap (2.7 eV). Following the band gap change as structural phase is changed, the 
onset photon energies of both the shift current susceptibility and the Glass coefficient are 
the lowest in the tetragonal phase, followed by that in the orthorhombic and rhombohedral 
phases, suggesting that the greater the lattice distortions, the higher the energy required to 
trigger a shift current response. It should be pointed out that the symmetric cubic phase 
is not able to exhibit shift current, although it has a lower band gap. Furthermore, for 
the same photon energies that are ~1.2 eV above the band gap of the tetragonal phase 
(£'°+1.2 eV, Eg is the LDA+t/ band gap of tetragonal KNbOs), all three phases exhibit 
comparable magnitude of both the shift current susceptibility and the Glass coefficient, with 
the response in the tetragonal phase slightly smaller than those in the other two phases. 

It is noteworthy that the band-edge (different Eg, and thus different photon energies) 
shift current response and Glass coefficient of the tetragonaf phase are much smafler than 
those of the other two phases. Figure shows the projected density of states (PDOS) onto 
the Nb d orbitals and the real-space wavefunction isosurfaces of the CBM at F A: point for the 
tetragonal and rhombohedral phases. Even though both the tetragonal and rhombohedral 
phases possess CBM composed of Nb d orbitals, they are subtly very different. The CBM 
of the tetragonal phase arises mainly from the Nb d^y orbital, with a 0.5 eV energy splitting 
between the d^y and dzx/zy orbitals. However, in the rhombohedral phase, there is almost 



1 0 1 
E-Ep(eV) 


FIG. 2. (Color online) The projected densities of states (PDOSs) onto the Nb d orbitals of the (a) 
tetragonal and (b) rhombohedral KNbOa. The inset shows the corresponding real-space wavefunc¬ 
tion distribntion aronnd the NbOe octahedra for the condnction band minimnm at the F point in 
the Brillouin zone. Upper: view along z] Lower: view along x. 
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FIG. 3. (Color online) The atomic structure representation of the (K,Ba)(Ni,Nb)05 solid solution 
with (a) layered and (b) rocksalt cation arrangements. Atoms are drawn with accepted ionic radii. 
In (a), both the apical and equatorial O vacancies are shown. 

TABLE II. The relative total energies (meV) of different O vacancy sites in tetragonal KBNNO. 
For the 1 x 1 x layered structure, the top apical site is in the K layer and the bottom apical site 
is in the Ba layer. 


Supercells Top-apical Bottom-apical Equatorial 


1x1x2 

+115 

+143 

0 

Rocksalt 

0 

0 

0 

1x1x3 

+94 

+124 

0 

1x1x4 

+77 

+103 

0 


no splitting between the dxy and d^x/zy orbitals, leading to a mixing of both dxy and dzx/zy 
orbitals in the CBM. This occurs because in the rhombohedral phase the Nb atom moves 
away from the Oe cage center along all three Cartesian directions, reducing the difference 
in the Nb-0 bonding strength in the different directions. In contrast to the dxy orbital, the 
dzx/zy orbitals extend the wavefunction along the shift current direction (z), facilitating the 
motion of the shift current carriers. Consequently, the shift current susceptibility and Glass 
coefficient are larger in the rhombohedral phase than in the tetragonal phase. 

B. Effect of O vacancy and cation arrangement on shift current in KBNNO 

In KBNNO, some of the original Nb^+ ions are substituted by Ni^+ ions (Ni^'^), with 
the charge compensated by the combination of A-site substitution of Ba^+ for K+ ions 


7 















(Bait) and O vacancies (V^o*). Generally, the O vacancies prefer to form adjacent to the 
Ni^"*" dopant, because the possible O vacancy as a donor is attracted by the acceptor. 
However, there are still three inequivalent O vacancy sites adjacent to Ni: the top apical 
site (with the NiJ((j,-V"o‘ local polarization parallel to the overall polarization), the equatorial 
site, and the bottom apical site (with the Nij((^-V"o* local polarization antiparallel to the 
overall polarization). We study extensively the stabilities of different O vacancy sites by 
using supercells with different compositions or cation arrangements, including both layered 
and rocksalt S-cation structures (Fig. |^. As shown in Table [TH in the 1 x 1 x A layered 
supercells, the equatorial O vacancy site is much more stable than both the top and bottom 
apical sites. This is because when the O vacancy is located at the equatorial site, a network of 
-Ni-Vo-Ni- forms. The attractive Coulomb interaction between the donor Vq and acceptor 
decreases the total energy, resulting in a more nearly uniform charge distribution than 
the other two cases. Also, the top apical site is slightly more favorable than the bottom 
apical site because of the repulsive interaction between the Vq and Baj^ donors. However, in 
the rocksalt arrangement, there is no preferred O vacancy site, as the charge environment is 
nearly isotropic for different orientations of the Ni-Vo complex. Even though the 
local polarization can be parallel or antiparallel to the overall polarization, it has minor effect 
on the preference of the O vacancy site. Basically, the rocksalt arrangement is about 100 
meV/atom less stable than the 1x1x2 layered arrangement with equatorial O vacancies, 
but more stable than that with apical O vacancies. Overall, we see a clear effect of the 
Coulomb interaction and charge compensation mechanism on determining the favorable O 
vacancy site. 

To study the effect of different oxygen vacancy sites and cation arrangements on shift 
current, we calculate the shift current of the 1x1x2 layered supercells both with equatorial 
and apical O vacancies, as well as the rocksalt B-cation arrangement. Although the layered 
structures with a high concentration of vacancies may be experimentally difficult to synthe¬ 
size under normal conditions, they serve as good examples to elucidate the stability of O 
vacancies and their effects on shift current. Figure]^ shows the calculated shift current sus¬ 
ceptibilities and Glass coefficients of these three different structures of the (K,Ba)(Ni,Nb )05 
solid solution. Several features are clear from the comparison of these results. First, the 
Glass coefficient of the layered arrangement with apical O vacancies (5xl0“^° m/V) is ap¬ 
proximately 12 times larger than that of the prototypical ferroelectric photovoltaic BiFeOs 




FIG. 4. (Color online) The largest tensor element of (a) shift current susceptibility {(Jzzz) and (b) 
Glass coefficient {Gzzz), and the (c) imaginary dielectric constant {ezz ), and (d) shift vector (Rz) 
of the (K,Ba)(Ni,Nb )05 solid solution with 1x1x2 layered and rocksalt cation arrangements. 
For the 1x1x2 layered arrangement, both the apical and equatorial O vacancy sites are included, 
while for the rocksalt arrangement the vacancy is at the apical site. 

(0.4x10“^° m/V) for the photon energies between its Eg and Eg+1.0 eV, but the required 
photon energies in KBNNO are much lower, at the LDA+U level of 1.0 eV. The low band 
gap of this KBNNO solid solution is further corroborated by its HSE06 band gap of 1.28 
eV (Table [^. Therefore, this cation arrangement would be a great bulk photovoltaic in a 
thick sample if the vacancy locations can be controlled. Compared to the Glass coefficient 
spectrum, there is not any peak in the shift current susceptibility spectrum or in the shift 
vector spectrum for photon energies near 1 eV. The Glass coefficient is large because the 
absorption coefficient at these energies is extremely small. These electronic transitions are 


9 




































FIG. 5. (Color online) The band structure of the 1x1x2 layered (K,Ba) (Ni,Nb )05 solid solution (a) 
with apical and (b) equatorial O vacancies. The inset shows the real-space wavefuction distribution 
for the corresponding conduction electronic states at ^4(0.5, 0.5, 0.5) k point in (a) and F k point 
in (b). The wavefunction plot in (a) is in view along x, while the left and right insets in (b) are in 
views along x and z, respectively. The K and Ba ions are not shown in the wavefunction plot. 


mainly from the O 2p orbital dominated valence band (VB) to the CB that is composed of 
Nb 4(i^2 and Ni 3 d ^2 orbitals as well as O 2p^. The transitions near the ^4(0.5, 0.5, 0.5) point 
of the Brillouin zone are not localized at the sublattice around Ni, but broadly distributed 
over the whole lattice [Fig. ga)]. 

Second, both the shift current susceptibility and Glass coefficient of the 1x1x2 layered 
KBNNO with apical O vacancies are much larger than those of layered KBNNO with equa¬ 
torial O vacancies for almost the whole spectral range. This is because both the imaginary 
dielectric constant and the shift vector of the layered apical vacancy structure are larger 
than those of the layered equatorial vacancy one. Figure shows the band structure of the 
1x1x2 layered KBNNO solid solutions with equatorial and apical O vacancies. Unlike 
the apical case, where the band-edge transitions are near the A point, in the equatorial 
structure they mainly occur near the F point. In the latter case, the CBM has a major 
contribution from the Nb Ad^y orbital [Fig. [^b)]. Unlike the ^^2 orbital, the d^y orbital is 
distributed within the xy plane that is perpendicular to the shift current direction (z), lead¬ 
ing to a much smaller onset Glass coefficient for the equatorial case. When the equatorial 
O vacancies are organized along the x direction, a chain of -Ni-Vo-Ni- formed. Since the 
remaining Ni 04 complex prefers a square planar symmetry in the perpendicular yz plane. 
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FIG. 6. (Color online) The projected densities of states (PDOSs) of the (K,Ba)(Ni,Nb )05 solid 
solution with 1x1x2 layered and rocksalt cation arrangements. The inset shows the real-space 
wavefunction distribution (view along x) for the corresponding electronic states as indicated by 
the red arrow. In both cases, the O vacancy is at the top apical site. In (a), the equatorial O atom 
nearest to the viewer is hidden in order to show the orbital character of the Ni Sd^y/zx orbitals. 


the overall lattice asymmetry along the z direction is significantly reduced. On the other 
hand, removing the O atom at apical site leads to a larger c lattice constant and an overall 
enhancement of the lattice asymmetry along the z direction. The combined effect of the 
different lattice asymmetries and the orbital character change dictates the change of the 
shift current magnitude. 

Also, comparison of the layered and rocksalt cation arrangements, both with the apical O 
vacancies, shows that the rocksalt cation arrangement exhibits a much larger shift current 
susceptibility, albeit with a higher onset photon energies, whereas its Glass coefficient is 
smaller at lower energies (<2.8 eV) but greater at higher energies (>2.8 eV) with respect to 
that of the layered arrangement (Fig. |^. Simultaneously, the imaginary dielectric constant, 
which is proportional to the transition rate, of the rocksalt cation arrangement is also smaller 
for photon energies smaller than 2.8 eV, but greater for the photon energies above 2.8 eV. 
Comparison of the shift vector shows that even though the rocksalt arrangement has a 
smaller shift vector for the whole spectral range, the difference in the magnitude is decreasing 
with increasing photon energies. Figure shows the PDOSs of both the 1x1x2 layered 
and rocksalt cation arrangements. Clearly, the top of the VB has Ni 2>dzy/zx orbitals in 
the layered arrangement, but the VB has Ni ‘idxy orbitals in the rocksalt arrangement, 
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leading to a preferred motion of the shift current carriers and a larger onset Glass coefficient 
in the layered arrangement for the electronic transitions at the band edge. However, as 
the absorbed photon energies increase and states below the VB becomes involved, Ni 3dz2 
orbitals become more important for the rocksalt cation arrangement while Ni Sd^y orbitals 
play a greater role for the layered arrangement. This leads to a steady enhancement of the 
shift vector magnitude in the rocksalt but not in the layered cation arrangement. 

Furthermore, the electronic states of both the VBs and CBs in the rocksalt cation ar¬ 
rangement are much more localized than those in the layered structure. This gives rise to a 
more sharply peaked DOS contributing to electronic transitions in a narrower range of pho¬ 
ton energies [Fig. |^c)], in agreement with the overall greater transition rate in the rocksalt 
cation arrangement. Combined, these two effects dictate the greater shift current suscepti¬ 
bility in the rocksalt cation arrangement. These features are not readily evident from the 
Glass coefficient, as the Glass coefficient inciudes the light attenuation effect represented 
by the absorption coefficient (imaginary dielectric constant) that relates to the strength of 
the electronic transitions. Consequently, the Glass coefficient of the rocksalt cation arrange¬ 
ment is only moderately larger than that of the layered cation arrangement and only for high 
photon energies. This more localized nature of the electronic states in the rocksalt cation 
arrangement can be ascribed to its structural properties. Compared to the layered arrange¬ 
ment, in the rocksalt cation arrangement the -Ni-Vo-Ni- network is interrupted by the Nb 
atoms, leading to more localized Ni 3d orbital states in the CB. This narrower bandwidth of 
the CB not only induces a larger band gap (Table and a higher shift current onset photon 
energy, but also an overall enhancement of the shift current magnitude. Therefore, we see 
that the resulting shift current is significantly affected by the orbital character of the elec¬ 
tronic transitions and the localization of these electronic states, which are in turn affected 
by the structural properties including the lattice asymmetry, the cation arrangement, and 
the location of the O vacancies. 

C. The shift current of KBNNO with a lower concentration of O vacancies 

The above solid solutions have a fairly high concentration of O vacancies, which could 
impede the motion of the photocurrent carriers, because they may behave as recombination 
centers. If so, it would be necessary to reduce the amount of O vacancies while preserving 
the beneficial effects of O vacancies in reducing the band gap and enhancing the visible- 
light absorption. In the following, we study the shift current of KBNNO with a lower 
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FIG. 7. (Color online) The Glass coefficient of the (l-x)KNb03-xBa(Nii/2Nbi/2)Oii/4 solid solu¬ 
tion with different compositions (x=2/3, 1/2, and 1/3). The left and right insets show the real-space 
wavefunction of the valence band maximum (VBM) of the l/3KNb03-2/3Ba(Nii/2Nbi/2)Oii/4 
{y/2 X ^/2 X 3) and l/2KNb03-l/2Ba(Nix/2Nbi/2)Oii/4 (v^ x x 4) solid solutions, respectively. 
The wavefunction in the right inset has an extended nature along the Cartesian z direction, while 
that in the left inset is only distributed within the xy plane. 

concentration of O vacancies. These KBNNO solid solutions have the compositions of (1- 
x)KNb 03 -xBa(Nii/ 2 Nbi/ 2 )Oii /4 (x=2/3, 1/2, and 1/3), with a vacancy concentration of 
5.6%, 4.2%, and 2.7%, respectively. The corresponding supercells are \/2 x \/2 x 3, \/2 x 
\/2 x 4, and 2x2x3, respectively. In each supercell, two Nb^+ ions are replaced with two 
Ni^+ ions, and the charge is compensated by the combination of an O vacancy adjacent to 
Ni and four K+ ions randomly substituted by the Ba^+ ions (BaJ^). 

Fig. HI shows the calculated Glass coefficient of different KBNNO compositions. All 
three solid solutions exhibit a maximum Glass coefficient with photon energies ~1.9 eV. 
Furthermore, the x=l/2 composition has the lowest onset photon energy, as its direct band 
gap is the smallest among all three KBNNO solid solutions. It also exhibits the largest Glass 
coefficient (3.8x10“^° m/V) for the photon energies below 3.0 eV, which is ten times larger 
than the maximnm Glass coefficient of BiFeOa (0.4x10“^° m/V). Similarly, the difference 
in the magnitnde of the Glass coefficient is attributed to the different orbital compositions 
of the contributing electronic states. For the a:=2/3 composition, the top of the VB is 
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predominantly composed of O 2 p orbitals combined with a slight contribution of the Ni and 
Nb dzy/zx orbitals. However, there is an extensive contribution of the Ni and Nb dz 2 orbitals 
to the top of the VB for x=l/2 (Fig. [^. This more extended wavefunction nature along the 
shift current direction allows for an easier motion of the shift current carriers, corresponding 
to the larger Glass coefficient for x=l/ 2 . 

Previous experiment has shown that the x= 0.1 KBNNO solid solution has the lowest 
band gap and exhibits the best photovoltaic performance |HS]- Therefore, in the following 
we choose the 2x2x3 supercell that gives our lowest calculated vacancy concentration of 
2.7% (x=l/3) to study the effect of different A- and B-cation arrangements on shift current. 
There are two Ni^+ ions, four Ba^"*" ions and one O vacancy in each supercell. First, the shift 
current is calculated with the four Ba^"*" cations distributed over 1,2, and 3 different layers 
(indicated as Hi, A 2 , and A 3 ) while the two Ni^"*" cations are kept at the body diagonal 
positions with respect to each other (Fig. [^. The A cation arrangement only has a slight 
effect on shift current, including both its magnitude and photon energies. This is because 
the ionic radii of the K+ (1.64 A) and Ba^"'' (1.61 A) ions are quite similar, and thereby the 
change of the distribution of the Ba^"'' ions has only a minor effect on the overall structure. 
Also, the valence state of the A cation has a very delocalized s orbital character, which also 
has a minor effect on the electronic structure. As a result, all three structures with different 
H-cation arrangements have nearly identical electronic properties, with an impurity state 
above the top of the VB arising mainly from the O 2 p and Ni 3d orbitals. However, when the 
H-cation arrangement is varied, there is a substantial difference in both the Glass coefficient 
magnitude and the photon energies that induce the largest Glass coefficient. Specifically, 
the B 3 cation arrangement has the lowest Glass coefficient, with the highest onset photon 
energy. The two Ni^"*" ions in this arrangement are aligned along the Gartesian .2 direction, 
with an O vacancy in between them, forming a Ni-Vo-Ni complex. Moreover, the Bi cation 
arrangement with two Ni^+ ions distributed along the body diagonal direction exhibits the 
largest Glass coefficient with a moderately large onset photon energy, whereas the Glass 
coefficient of the B 2 cation arrangement is the second largest, but its onset photon energy 
is the lowest. 

We choose the Bi and B 3 cation arrangements as typical examples to study the underlying 
origin of the difference in their shift current responses. PDOS analysis shows that the Bi 
and B 3 cation arrangements have very different electronic structure properties. The Bi 
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FIG. 8 . (Color online) The atomic structure representation of the 2 / 3 KNb 03 - 
l/ 3 Ba(Nii/ 2 Nbi/ 2 )Oii /4 (2x2x3) solid solution with (a) two Ni^+ ions at the body diagonal 
positions {Bi) (b) two Ni^+ ions along the Cartesian 2 ; direction (S 3 ). The atoms are shown with 
their accepted atomic radii. The K+ and Ba^+ are omitted for clarity. The Glass coefficients of the 
2 / 3 KNb 03 —l/ 3 Ba(Nii/ 2 Nbi/ 2 )Oii /4 solid solutions with different (c) A-site and (d) S-site cation 
arrangements. The S-cation arrangement remains the same (with two Ni^^ ions distributed along 
the body diagonal direction) when the A-cation arrangement is varied in (c), whereas the A-cation 
arrangement remains the same when the S-cation arrangement is varied in (d). “Ai”, “A 2 ”, and 
“A 3 ” correspond to configurations with the four Ba^+ cations distributed within 1, 2, and 3 layers, 
while “Si”, “S2”, and “S3” correspond to the configurations that the two Ni^+ cations are dis¬ 
tributed along the body-diagonal, face-diagonal, and the Cartesian 2 : directions. “A 2 ” and “Si” 
are the same structure. 


arrangement has a gap state just above the top of the VB, leading to a smaller band gap 
and lower shift current onset photon energies with respect to those in the S 3 arrangement. 
Furthermore, there is a predominant Nb and Ni dz ‘2 orbital character in the Si arrangement. 
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FIG. 9. (Color online) The densities of states (DOSs) of the 2 / 3 KNb 03 —l/3Ba(Nii/2Nbi/2)Oii/4 
solid solutions with (a) the two Ni^+ cations aligned along the body-diagonal direction (Bi) and (b) 
the two Ni^+ cations aligned along the Cartesian z direction (S 3 ). The inset shows the real-space 
wavefunction distribution for the top of the VB state, which shows a great contribution of Ni 3d^2 
(also some Nb 4d^2) orbital character. 


but in the S 3 arrangement the O 2 p orbitals make the major contribution (Fig. |^. The 
extended wavefunction nature along the z direction is beneficial for the shift current 
response, resulting in a larger Glass coefficient in the Si arrangement. This difference in the 
electronic properties is ascribed to their different structural properties. Basically, there are 
two Ni^+ ions and one O vacancy in each supercell. In the Si arrangement, these two Ni^+ 
ions have different crystal environments: one with five adjacent O atoms (NiOs), the other 
with six (NiOe). However, in the S 3 cation arrangement, the two Ni^+ ions share the same O 
vacancy, corresponding to a NiOs environment for both Ni^^ ions. For the octahedral SOe 
complex, the d orbitals of the S cation split into triply-degenerate t 2 g and doubly-degenerate 
eg states. Because the Nb^+ ion is strongly ferroelectric, its off-center displacement leads 
to a concurrent change of the Ni-0 distance along the ^ direction for the NiOe complex. 
This induces additional splittings between the and dx 2 _y 2 orbitals of the Cg state and 


between the d^,, and d. 


orbitals of the t 2 g state (Fig. 10). This change of the structural 


asymmetry is only moderate. However, for the Ni 05 complex, the removal of one O atom 
at the apical site provides space for the O atom at the opposite site to move away from 
the central Ni atom. Correspondingly, the Ni-0 distance along the Cartesian z direction 
is much larger than that in the xy plane, giving rise to a near-square-planar symmetry, as 
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FIG. 10. (Color online) The schematic representation of crystal field splitting of the Ni^"*" 3d 
orbitals nnder different crystalline environments. The Ni-0 distance along the Cartesian z direction 
is longer than that in the xy plane. 

shown in Fig. [T^ Compared to the distorted NiOe environment, the splitting between the 
dz 2 and dx 2 _y 2 orbitals is much larger for the NiOs complex. Therefore, the energy of the ^^2 
orbital in the NiOs complex is much lower than that in the NiOe complex. As a result, the 
higher energy of the d^^ orbital in the NiOe complex induces a gap state for the Bi cation 
arrangement. The d^^ nature of the gap state generates a larger shift current response, as 
indicated by the Glass coefficient (Fig. [^. 

The calculated maximum shift current susceptibility of the 2x2x3 solid solution is 
1.5x10“^ Acm“^/Wcm“^, which is six times as large as the experimental observation for 
KBNNO (the samples are not completely poled) j35], and is also comparable to that of 
BiFeOa. Experiment has found that the short-circuit photocurrent of KBNNO is 0.1 /rA/cm^ 
under 4 mW/cm^ illumination with above-band-gap light, corresponding to a current den¬ 
sity of 0.25x10“^ Acm“^/Wcm“^, while for BiFeOs, the current yield is 4 jiK/cw? under 10 
mW/cm^ illumination with green light, for a current density of 4x 10 ^ Acm ^/Wcm ^ [21 
[35] . If we include the light attenuation effect in a thick sample of KBNNO, the photocurrent 
evolution is 0.4 mA/cm^ for a 100-nm-thick sample under the illumination of a 1000 W/m^ 
solar simulator, as estimated by Jz=GzzZ x h/d, where d and Iq are the sample thickness 
and light intensity. 
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FIG. 11. (Color online) (a) The Glass coefficient of the 2 / 3 KNb 03 -l/ 3 Ba(Nii/ 2 Nbi/ 2 )Oii /4 solid 
solution with different in-plane compressive strains, (b) The projected density of states (PDOS) 
of 2 / 3 KNb 03 -l/ 3 Ba(Nii/ 2 Nb 4 / 2 )Oii /4 solid solution under 2% in-plane compressive strain. The 
inset in (b) shows the real-space wavefuction distribution (side view) for the top of the VB state. 
The solid solution is with the “Bi” arrangement; the two Ni atoms are distributed along the body 
diagonal positions and the O vacancy is at the apical site of one Ni atom. 


D. The effect of strain on shift current 

It has been shown that strain can substantially affect the octahedral cage distortions, 
rotate the polarization, and change the band gap of perovskite oxides [22l [521, [53]. Each of 
these can have a significant effect on the material’s shift current response. In order to study 
the effect of strain on shift current, we apply in-plane compressive strains to the 2 / 3 KNb 03 - 
l/ 3 Ba(Nii/ 2 Nbi/ 2 )Oii /4 solid solution, in which the two Ni atoms are distributed along the 
body diagonal positions (Bi). In this configuration, the O vacancy is adjacent to only one 
Ni atom, leaving an O 5 environment around this Ni, but an Oe environment to the other Ni. 
Figure shows the calculated Glass coefficient as a function of the in-plane compressive 
strain. It is evident that the magnitude of the Glass coefficient decreases steadily with 
increasing in-plane compressive strain. Goncurrently, the shift current onset photon energy 
also increases, suggesting a bigger band gap with enhanced strains. This is contrary to the 
naive expectation that the magnitude of the Glass coefficient will increase with strain because 
of the polarization rotation towards the [ 001 ] direction and the overall enhancement of the 
structural asymmetry along the z direction when applying in-plane compressive strains. 

However, this change can also be rationalized by the electronic structure and wavefunction 
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analysis. PDOS analysis shows that the impurity state that is originally above the top of the 
VB shifts downwards and finally merges into it as the in-plane compressive strain increases 
[Figs. i a) and [IT|[b)]. Concurrently, the VB edge becomes predominantly composed of O 
2p orbitals [Fig. in comparison to the significant contribution of the ^^2 orbital for 

the configuration without strain. This not only gives rise to a larger band gap, but also 
affects the motion of the shift current carriers, leading to the observed reduction of the 
Glass coefircient. Further examination of the structure shows that the difference between 
the Ni-0 distance in plane and that along the .2 direction is substantially enhanced when 
the in-plane compressive strain increases. This results in a larger splitting between the ^^2 
and dx 2 -y 2 orbitals for the Ni with an Oe environment, becoming similar to the Ni atom in 


an O 5 environment (Fig. 10). This shifts downwards the original orbital dominated gap 
state, giving rise to both a larger band gap and a smaller Glass coefficient. 

IV. CONCLUSIONS 


In summary, we study the bulk photovoltaic effect of the prototypical KNbOa and visible- 
light ferroelectric photovoltaic (K,Ba)(Ni,Nb) 03_5 from first principles. The effect of lattice 
distortion, oxygen vacancies, cation arrangement, composition, and strain on shift current 
are systematically studied. We find that the maximum shift current response (with UV 
absorption) of the orthorhombic KNbOs is more than twice that of BiFeOs, although KNbOs 
has a wider band gap. Furthermore, the band-edge shift current response of the tetragonal 
KNbOs is smaller than that of its rhombohedral counterpart. This occurs because a more 
isotropic lattice distortion in the rhombohedral phase reduces the splitting between the d^y 
and dzx/zy orbitals, leading to a larger electronic contribution of the 2 :-direction-extended 
state to the shift current. 

In (K,Ba)(Ni,Nb) 03 _a, the charge compensation affects significantly the favorable loca¬ 
tion of oxygen vacancies, with more stable O vacancy location based on more effective charge 
compensation. The layered arrangement with apical O vacancies exhibits a much larger shift 
current response than that with equatorial O vacancies, as the lattice asymmetry is enhanced 
in the former, but reduced in the latter case. Compared to the layered arrangement, the 
rocksalt arrangement has more localized electronic states, giving rise to a larger density of 
electronic transitions within a narrower energy range. Combined with the effect of the more 
extended d^^. orbital, this gives rise to a larger shift current susceptibility. The effect of 
composition and A-cation arrangement on shift current is moderate, whereas the B-cation 
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arrangement substantially alfects both the electronic structure and the shift current, which 
can be rationalized by crystal field theory analysis. With an Oe environment around some 
of the Ni, the band gap is lowered, and the final shift current yield is enhanced. The ef¬ 
fect of strain on shift current is indirect, through the change of the wavefunction nature 
of the contributing electronic states, which can be used to engineer the shift current in a 
predictable fashion. It is noteworthy that the shift current response of (K,Ba)(Ni,Nb) 03_5 
is comparable to that of BiFeOs, but at a much lower photon energy. More importantly, 
the order of magnitude enhancement of the shift current response can be gained by simply 
layering the (K,Ba)(Ni,Nb )05 solid solution, with its Glass coefficient reaching 12 times that 
of BiFeOg. 

Finally, we have demonstrated extensively that the shift current is dictated by the wave- 
fuction nature of the contributing electronic orbitals, which are in turn unanimously related 
to their structural properties. Therefore, we have built up a bridge between materials’ struc¬ 
tural properties and their photovoltaic performances, and provided a pathway for analyzing 
and elucidating the connections among these different physical properties. The elucidated 
connection between materials’ structure, electronic structure, and shift current is useful for 
the design of bulk photovoltaic materials and understanding their PV mechanisms. 
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